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We study complex eigenvalues of large N x N symmetric random matrices of 
the form TC — H — iT, where both H and F are real symmetric, H is random 
Gaussian and f is such that iVTrf \ ~ TrH? when N -> oo. When f > the model 
can be used to describe the universal statistics of S-matrix poles (resonances) in 
the complex energy plane. We derive the ensuing distribution of the resonance 
widths which generalizes the well-known \ 2 distribution to the case of overlapping 
resonances. We also consider a different class of "almost real" matrices when T is 
random and uncorrelated with H . 
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As is well-known, S-matrix poles or resonances can be looked at as complex eigenvalues of an 
' effective non-Hermitian Hamiltonian H emerging after the so-called complex rotation was applied 
to the original Hermitian Hamiltonian of the system. The method JjJ allows one to extract reso- 
nance positions and widths T^ directly, without expensive evaluation of the energy-dependent 
S— matrix elements. The search of resonances typically amounts to diagonalizing large complex 
symmetric matrices representing the effective Hamilonian in a suitable basis on a finite grid. Us- 
ing this method, patterns consisting of hundreds of individual resonances were extracted recently 
for realistic models in atomic and molecular physics ||||. Frequently, for high enough excitation 
energies the S-matrix poles are placed irregularly in the complex plane, forming a structure of 
a "chaotic jumble" ||, the situation calling for a statistical description. As another important 
development it is appropriate to mention a powerful numerical algorithm that was proposed to 
extract resonance positions and widths with a high accuracy from the experimentally accessible 
time correlation function of finite duration All these facts make the problem of statistical 
description of resonances to be an important and opportune task. 

At the same time, the statistics of highly excited bound states of closed quantum systems of 
quite different microscopic nature is known to be system-independent (universal), provided the 
corresponding classical counterparts demonstrate a well-developed chaotic motion J5j. Moreover, 
the spectral correlations turn out to be exactly the same as provided by the theory of large random 
matrices on the scale determined by a typical separation A between neighbouring eigenvalues. 
Microscopic justifications of the use of random matrices for describing the universal properties of 
quantum chaotic systems have been provided by several groups recently, based both on traditional 
semiclassical periodic orbit expansions and on advanced field-theoretical methods ft is 

natural to try to develop a description of universal statistical properties of resonances using the 
same ideas. 

The methods to adjust random matrix description to the case of open chaotic systems are well 
known since the pioneering paper [fio| and described in much detail in The starting point of 
this approach is the representation of the scattering matrix S in terms of an effective non-Hermitian 
Hamiltonian TLeff = H — iT. Here the Hermitian Hamiltonian H describes the closed counterpart 
of the open system whereas the anti-Hermitian part iT arises due to coupling to open scattering 
channels. It has to be chosen in the form ensuring the unitarity of the scattering matrix. As a 
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result all eigenvalues of f turn out to be non-negative (as is also required by causality), the number 
of strictly positive eigenvalues being just the number M of open channels. 

It is natural to expect, that statistical properties of resonances are inherited from the correspond- 
ing random matrix universality of levels of closed systems. Remembering the discussion above, 
we expect them to be universal on the energy scales in the complex plane comparable with the 
mean level spacing for the closed system A. In contrast, properties on a much larger scale must 
be highly system-specific. 

A natural way to incorporate the random matrix description of the quantum chaotic system is 
to replace H by a large N x N random matrix of appropriate symmetry. Namely, chaotic systems 
with preserved time-reversal invariance (TRI) should be described by matrices Hij which are real 
symmetric. Such matrices form the Gaussian Orthogonal Ensemble (GOE), whereas for systems 
with broken TRI one uses complex Hermitian matrices from the Gaussian Unitary Ensemble (GUE) 
||. As to the matrix f it is determined by the coupling of the chaotic system to open channels 
and can be considered as a fixed given one of the same dimension. 

In general, the two matrices H and T do not commmute, the fact making the analysis of complex 
eigenvalues of H. to be a rather non-trivial problem. In the simplest case one can treat the anti- 
Hermitian part as a small perturbation to the Hermitian one. This is justified for the regime of 
isolated resonances, i.e. when typical resonance widths are much smaller than the mean separation 
A between the positions of neighbouring resonances. Assuming for simplicity that all scattering 
channels (hence, the eigenvalues of T) are statistically equivalent, one arrives at the so-called x 2 
distribution of the resonance widths. It was introduced originally many years ago by Porter and 
Thomas and since then experimentally observed in many physical systems, see some references 
in P2| . The corresponding perturbative expression is known also for the case of non-equivalent 
scattering channels JT|| . 

Until quite recently relatively little was known about resonance statistics non-perturbatively. 
Among the known facts were the joint probability density of all resonances (however, not its 
moments) for the systems with only one open channel M = 1 jl^l , and the mean density of 
complex eigenvalues for very many open channels M ~ -/V 1 Jig] , The most essential progress 
was achieved recently by employing a powerful field-theoretical method known as Efetov's 

supersymmetric non-linear a— model [ p~7| . It turns out to be possible to find explicitly the mean 
density of the resonances for systems with broken time-reversal invariance and to put forward a well- 
grounded conjecture on the higher correlation functions in the complex plane fl8|| . Unfortunately, 
the most interesting and practically important case of the systems with preserved time-reversal 
invariance proved to be less accessible analytically (for numerics see JlS|] ) and no results of such 
generality were obtained for such systems so far. 

The main goal of the present paper is to fill in this gap partly by presenting the exact non- 
pertubative distribution of the resonance widths for quantum chaotic systems with preserved time- 
reversal invariance. 

It is necessary to menion that random matrices (more generally, random linear operators) with 
complex eigenvalues emerged recently in many other physical contexts. A few recent examples 
include the description of thermal motion of an isolated vortex in disordered type-II superconduc- 
tors with columnar defects , the problem of a classical diffusing particle advected by a stationary 
random velocity field and attempts to understand the universal features of chiral symmetry break- 
ing in Quantum Chromodynamics (a more detailed discussion and the relevant references can be 
found in the papers mentioned below). At the same time, our knowledge of statistical properties 
of random non-self adjoint matrices is quite scarce and incomplete. This fact recently stimulated 
efforts of different groups to improve our understanding in this direction |2C|-|2"4f| . In particular, the 
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existence of a nontrivial regime of weak non-Hermiticity was recognized in the work p2fl , see more 
detailed discussion in H] . The guiding idea to realize the existence of such a regime comes from 
the experience with resonances |]l6| . Guided by that example one guesses that a new regime occurs 
when the imaginary part of typical eigenvalues is comparable with the mean separation between 
neighboring eigenvalues along the real axis. 

One can use again the supersymmetry approach to obtain the mean density of complex eigen- 
values in the regime of weak non-Hermiticity for matrices with independent elements |2^ , |24j| . The 
density turned out to be described by a formula containing as two opposite limiting cases both the 
Wigner semicircular density of real eigenvalues of the Hermitian matrices and the uniform density 
of complex eigenvalues discovered for usual non- Hermitian random matrices already by Ginibre 
p5| , in much details addressed by Girko (2(| and other authors [^7| . In his insightful paper Efetov 
p8f managed to derive the density of complex eigenvalues for a related, but different set of (almost- 
symmetric) real random matrices. His result was generalized in |Tq] . Actually, the non-Hermitian 
matrices considered in p^ ] and |28j are just two limiting cases of a general three-parameter family 
of "weakly non- Hermitian" ensembles |^4| with independent entries. Apart from these two cases, 
the mentioned three-parameter family contains one more important particular case corresponding 
to an ensemble of complex symmetric random matrices. Such matrices are natural to call the 
"almost-real" ones. The second goal of the present paper is to provide the eigenvalue density for 
this particular ensemble. 

As a starting point of our investigation we use the representations of the mean density p(X, Y) 
of eigenvalues Z = X + iY with a given imaginary part Y derived by us earlier |24|,^{|. Let 
us first consider the case relevant to the chaotic scattering. Due to rotational invariance of the 
GOE matrices the result can be dependent only on the set of positive eigenvalues 71,72, ...,7m of 
the matrix T. Correspondingly, the density px{y) = p(X,Y)A 2 (X)/n of scaled resonance widths 
y = (measured in units of the local mean level spacing A(A) = l/{Nv{X)) of the closed 
system, with v{X) = 7r _1 (l - A 2 /4) 1 / 2 being the semicircular density of real eigenvalues) for those 
resonances whose positions are within the narrow window around the point X of the spectrum is 
given by: 

Mr.-,. 

(px(y)) = ^ / dMQ)Str(4 F )Q)str(a T Q) Sdet" 1 / 4 1- ^~{Q,^} exp ^Str (cr T Q) 

' J o=l _ 9a _ 



(1) 



We introduced here the parameters g a G [1, 00) characterizing the coupling of a system to open 
channels which are expressed in terms of the eigenvalues 7 a as g a = 2 ttv(X) + la 1 )- The 
integration above goes over the set of 8 x 8 supermatrices Q satisfying the constraint Q 2 = — 1 
and some symmetry requirements, which can be conveniently resolved if one parametrizes the 
manifold as: Q = —iZAZ^ 1 , with A = diag(l, — 1), each 1 representing 4x4 unity matrix. The 
symbols Str, Salet stand for the graded trace and the graded determinant, correspondingly, and 
|<5,<5y j = Qa T + <j t Q stands for the anticommutator. Other 8x8 supermatrices entering this 
expressions are as follows: 
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and T3; ,fg are 4x4 diagonal supermatrices: f 3 = diag{f,f}; fg = diag{02,f}. with f = 
diag(l, -1). 

For the case of the matrix T being real symmetric with entries independently fluctuating around 
zero and normalized in such a way that limjv^oo Trr 2 < 00 (the condition of being "almost 
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symmetric") the corresponding expression for the mean density of eigenvalues is obtained from 
Eq.(Q) by replacing the product of superdeterminants in the integrand with the exponent: 

„2 



exp j-^j-Str (a T Q} 



where v 2 = (ttv) 2 limjv^oo Trf |, see the derivation in [ p4| . 

The starting formulas are actually valid both for real symmetric and complex Hermitian matrices 
H. To extract the explicit form of the distribution function one still has to perform the integration 
over the manifold of the supermatrices Q. Still, it is rather difficult calculation due to a cumber- 
some structure of that manifold, and the success is mainly determined by finding an appropriate 
parametrization. For the simplest case of matrices H being complex Hermitian the supermatrices 
Q become block-diagonal due to additional constraints and the calculation, although cumbersome, 
can be still done successfully [^6| , p]j2^1 in the "standard" Efetov parametrization jl7]]. For the 
case of real symmetric matrices with resulting full 8x8 structure of Q this is no longer the case. 
Fortunately, we are able to parametrize the matrices Z as Z1Z2 in such a way, that Z\ commutes 
with <t t , thus dropping out from the integrand. We would like to stress that although we enjoyed 



some insights from the recent paper on real asymmetric matrices 28 , our actual parametrization 
is quite different from that used in |28[ . A derivation of the explicit form of such a parametrization 
and especially the integration measure turn out to be very lengthy and will be presented in details 
elsewhere. Here we would like to present only the result of integrating out the Grassmannian 
degrees of freedom. 

where we use the convention s = r for the case of resonance statistics and s — R for that case 
when the antihcrmitian part F is a random matrix. Correspondingly, we defined 

e-scn \ii ( TTm — for resonances, s = r 

(1 _ A 2)l/2 

exp (v 2 X 2 ) for a random matrix, s = R 

The following important comment is appropriate here. The integrals over Ai, A2 as they stand 
in Eq.(Q) are not well-defined when taken literally along the real axis due to a singularity at 
the points Ai = A = 0; A2 = A = 0. This difficulty can be traced back to a singularity of 
the transformation arising from diagonalization of some supermatrices when coming to a new 
parametrization mentioned above. A more accurate consideration shows that the integration over 
Ai, A2 should be understood in the sense of principal value. To be precise, the correct expression 
is equal to the half-sum of two integrals with integration contours encircling the singular point 
Ai = — iX; X 2 — —iX from above and from below. 

The expression Eq.(||) is exact, but still is not very convenient for evaluation. Fortunately, for the 
physically motivated case of resonances one can deform the contour in such a way, that it goes along 
the cut A1.2 = —ip, V £ [1> 00 ) f° r negative values of y and along a similar cut A1.2 = ip, p G [1, 00) 
for y > 0. At the same time, all the singularities Ai = —ig a , a = 1, ...,M belong to the cut in 
the lower half-plane. This analytic structure of the integrand ensures vanishing of the density for 
y > as expected, whereas for y < one obtains: 

1 d 2 f 1 2 f°° dpi f pi d P2 (pi - p 2 )e+y^ 1+ P 2 - 2 ^ 



M , ., (a) 



x xi(p2)x2(pi) TT ?i ~i\ — ~ — n~' 

„ = 1 (\9a -Pl\\9a -P2\) 1/2 
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where the functions Xi.2{p) are defined as follows: 

1 for 94k <P< 94k+i f 1 for ,94fc+i < P < 9ik+2 

Xi{p) = { -1 for g ik+2 < P < 94k+3 Xa(p) = I -1 forg 4fc+3 < p < g 4k +4 
otherwise I otherwise 

where k = 0, 1, [M/4] and we used a convention: 1 = go < gi < ... < gu < 9m+i = oo. 

The formulas Eqs.(^|) provide us with the most general explicit analytical expression for the den- 
sity of complex resonances for a chaotic quantum system with preserved time-reversal invariance. 
As such they constitute the main result of the present paper. 

We see, that in general these formulas are still quite cumbersome. One can infer from them 
various known limiting cases considered earlier. For example, we can consider the case of weak 
coupling to continuum when all g a 1. By noticing that the integration over p2 in Eq.(|^) is 
dominated in this case by p2 ~ 1, whereas p\ > gi » 1,X,P2 and taking into account also that 
typical values of y are of the order of g^ 1 << 1 one can evaluate the integrals over p2, A to the 
leading order. The result is: 

px(y < 0) = -Il5a /2 r rr i ~ n/2 e ^2(pi) (4) 

which coincides with the distribution derived recently by Alhassid and Lewenkopf |l3[ ] by a per- 
turbative treatment, up to alternating signs in the factor xz(p) which is an apparent misprint in 
their paper. 

Relying on earlier experience with resonance statistics in systems with broken time-reversal 
invariance |16| , |rl| ], one might expect that the distribution simplifies drastically for the case of 
statistically equivalent channels: g = gi =52 = ••■ = 9m- Surprisingly, this seems to be not 
the case for the present model. Actually, we find it difficult to perform such a limit explicitly 
for arbitrary number of open channels, and even for M = 4 the result turns out to be quite 
cumbersome. Below we present the formula for one and two open channels M = 1,2: 

Px(y < 0) = J d\{l - A 2 ) e - 2 ^F M (A, y) (5) 



(Pi -P2) 



where for M — 1 we found 

r°° e ypi rs e yp2 

Fl = {g ~ X) ] g dPl (pi-lVmX-p 1 n Pl -g)^J 1 dP \pl-l)^\X-P2Y(9-P2)^ 

and the case M = 2 is the simplest one: 

7T f 9 1 e y(P2+g) 

F " - ( 5 2 - 1)1/2 ] 1 P \pl-iy/2 (A-p 2 )2 

The figures Fig.l and and Fig. 2 show that the resulting expression favourably agree with the 
results of direct numericaldiagonalization of the corresponding complex matrices. 

Despite the fact, that the general resonance widths distribution for the case of preserved TRS 
turns out to be much more complicated than for the case of broken TRS, both expressions share 
many important features. For example, evaluating the first moment of the distribution Eq.(|j|) for 
M = 2 exactly, one arrives at the following expression for the mean resonance widths: 

A 7T .9+1 
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This formula is well known in nuclear physics as Moldauer-Simonius relation and follows from the 
unitarity of the scattering matrix ]30| . Actually, a simlar relation must be satisfied for an arbitrary 
number M of open channels. It was indeed found to be the case for systems with broken TRS 
. Unfortunately, it is quite difficult to verify the Moldauer-Simonius relation directly from the 
general distribution Eq.(||). 

The logarithmic divergency at g a = 1 signals on the powerlaw tail behaviour px(\y\ 3> 1) oc y~ 2 
typical for this (a so-called "perfect coupling") case. It is instructive to perform the limit 
directly in Eq.(|) for M = 2 and find: 

p x ( y< Q) = -L (2y + l + (2y-l)e 4 y) (7) 

where such a behaviour is self-evident. 

Actually, it is interesting to compare this expression with the M = 1 case for systems with 
broken TRS 

Px(y<0) = ^(l + (2y-l)e 2 y) (8) 

The reason is that in the regime of isolated resonances g ^> 1 the M-channel system with broken 
TRS possesses exactly the same widths distribution as the 2M channel system with unbroken TRS. 
A similar correspondence also holds in the limit of many open channels M ^> 1. Moreover, before 
the exact formula Eq. (R) was available, an attempt was made J31J to compare M = 1 formula Eq. (||) 
valid for systems with broken TRS with numerically found resonance widths in M = 2 chaotic 
TRS-preserving microwave cavity. The correspondence was reported to be quite satisfactory. 

We see, that both formulas indeed give exactly the same behaviour for the resonances whose 
widths exceed the typical separation between them: px(|y| 3> 1) ~ l/2y 2 . However, in the region 
of very narrow resonances the behaviour is slightly different: px(y — > 0) = 1 for broken and p(y — > 
0) = 4/3 for unbroken TRS. This fact shows the limitations of the mentioned correspondence. 

As to the random-perturbation case (s = R in cq.(|^), the presence of the Gaussian factors 
in the integrand of Eq.(^|) prevents one from deforming the contour to the complex plane. At 
the moment, we are unable to simplify the expression further. However, one can still infer that 
in the perturbative case d<1 the function px(y) reduces to a simple Gaussian, whereas in the 
opposite limiting case b>1 one indeed arrives at the "elliptic law" well-known to hold for strongly 
non-Hermitian random matrices . 

The authors are much obliged to B.A.Khoruzhenko for useful communications. The work was 
supported by SFB 237 "Unordnung und grosse Fluktuationen" and in part by RFBR Grant No. 
96-02-18037-a. 
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FIG. 1. Comparison between the predictions of our formula Eq.(g) for one channel case M=l at "perfect 
coupling" g=l and the results of direct numerical diagonalization of K = 10000 complex matrices of the 
size 400 x 400. 
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FIG. 2. Same as Fig.l for M=l, g=3. 



